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Abstract 

This work reports the structural instability at T=0 °K of the U 2 M 0 compound in the structure Cl U under the distortion 
related to the Cee elastic constant. The electronic properties of U 2 M 0 such density of states (DOS), bands and Fermi 
surface (FS) are studied to understand the source of the instability. The Cllt structure can be interpreted as formed 
by parallel linear chains along the z-directions each one composed by successive U-Mo-U blocks. The hybridization 
due to electronic interactions inside the U-Mo-U blocks is slightly modified under the distortion. The change 
in distance among chains modifies the U-U interaction and produces the split of f-states. The distorted structure is 
stabilized by the energy lowering of the hybridized states, mainly between d-Mo and f-U states, together with the 
f-band split. Consequently, an induced Peierls distortion is produced in U 2 M 0 due to the distortion. In addition, 
the results of this work remarks that the structure of the ground state of the U 2 M 0 compound is other than the assumed 
Cl lb structure. 


1. INTRODUCTION 

The bcc U-Mo solid solution has technological relevance because is a promising candidate for developing high- 
density uranium fuel, to be used in research reactors with low enrichment uranium (LEU, 235U < 20 at%) [1]. It 
could be also a potential fuel for transuranic-burning advanced reactors [2] . 

The development of high-density uranium fuels is one of the major challenges found in Material and Test Reactors 
research. Of particular interest is the U-Mo solid solution in the y-phase dispersed in an A1 matrix or in a monolithic 
form. Both of them generated high expectation several years ago due to its apparent acceptable irradiation behaviour, 
low to moderate fuel/matrix interaction, and stable fission gas bubble growth for moderate neutron fiuxes. However, 
unexpected failures like pillowing and large porosities in LEU UMo dispersion plates in high neutron flux irradiation 
experiments have been reported. Since these failures were found in 2004 [3-5], much effort has been devoted to the 
search for their root causes and, mainly, to find a solution. It is now generally accepted that the properties of the 
interaction layer formed between the (U, Mo) fuel and the A1 matrix are important factors in the failures. Conse¬ 
quently, the research has been focused in the development of diffusion barriers in the interface by adding additives 
to the fuel and/or to the matrix, mainly Si additions to the A1 matrix. However, from the experimental point of view, 
it appears nowadays that the addition of Si in itself may be insufficient to meet the requirements of relatively high 
power at high burnup [6] . In addition, theoretical results using BFS methods for alloys showed that the formation of 
efficient diffusion barriers is strongly dependent on the Mo concentration near the interface [7]. In fact, the analysis 
shows that Si diffuses deeply in the U-Mo bulk for low levels of Mo in areas close to the interface. As the minimum 
necessary amount of Mo near the interface can not be controlled in a solid solution, an effective diffusion barrier can 
only be achieved locally in areas with enough Mo content near the interface. This theoretical result suggests that the 
possibility of obtaining a reliable and licensable fuel based on (Al,Si)/(U,Mo) could be a endless line of research [8]. 
The theoretical modelling also provides alternatives to the Si addition such as the C coating to produce an absolute 
diffusion barrier [9] . In any case, it is important to remark that the still unsuccessful development of (U, Mo) fuels 
leads to the need for other alternatives. Therefore experimental and theoretical research is needed to develop a reliable 
fuel. Besides the modelling effort using the BFS method, there is no other theoretical works known by the authors of 
this work aimed to describe the atomic behaviour in the interface. A few theoretical efforts using ab initio methods 
have been made [10-12] to understand the main thermodynamic properties of the U-Mo alloy and the stability of the 
U 2 M 0 compound in the Cl 1^ (MoSi 2 prototype) structure, which is the only stable compound observed in the phase 
diagram. The importance of ab initio works is the possibility to describe the main physical properties of a system 
with no experimental information. However, the ground state must be properly characterized, mainly, if the cluster 
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expansion technique is applied to obtain the correct behaviour of the disordered U-Mo alloy [13-15]. The work of 
Alonso et al. [11] studied the formation energy of the U-Mo solid solution in the bcc phase using the ah initio cluster 
expansion technique with a set of effective cluster interactions (ECI). The ECI were computed using the direct in¬ 
version method from the total energies of 16 bcc-based ordered structures calculated by the Eull-Potential Linearised 
Augmented Plane Wave method (PPLAPW)[16]. These authors calculated the formation energy of the disordered y 
U-Mo solid solutions as well as the U 2 M 0 (CIU) compound. The calculated energies of formation predict only one 
stable compound (U 2 M 0 ) and an asymmetry that stabilizes the disordered alloy on the uranium rich side at 0°K. In 
addition, temperatures of calculated phase equilibrium using CVM were excessively high. The calculated disordering 
temperature of the U 2 M 0 (CIU) compound is ~ 2000°C, higher than the experimental value of ~ 600°C [17]. This 
failure indicates that more clusters should be included in the cluster expansion and a bigger maximum cluster should 
be selected in the CVM method to study phase equilibria in U-Mo system. 

A more precise calculation of the heat of formation, the bcc ground state and phase diagram is presented by 
Landa et al. [10]. These authors applied the Green’s-function technique based on the improved screened Ko- 
rringa-Kohn-Rostoker method combined with the coherent potential approximation (CPA) to treat compositional 
disorder. In addition, the authors used the EPLMTO methods to compute the heat of formation of compounds and a 
16-atom supercell to treat the compositional disorder within this formalism. Calculated heats of formation are com¬ 
pared with CALPHAD assessments. The decomposition curves for y-based U-Mo solid solutions are derived from 
Ising-type Monte Carlo simulations. These authors discussed how the heat of formation correlates with charge transfer 
between the alloy components and found that the stabilization of the Cl U-phase in the U-Mo system is explained in 
terms of the electronic DOS change due to ordering of the U 2 M 0 alloy, i.e. a significant drop of the DOS in the vicinity 
of the Eermi level in the case of the ordered Cl U compound causes a decrease of the total energy. Due to significant 
size mismatch between U and Mo atoms, these authors emphasize the importance of including the compositional and 
elastic dependence of ECI to produce lower temperature of decomposition of the y U-Mo solid solutions. 

Remarkably, at the present time, only a few experimental data have been published related mainly to thermal and 
mechanical properties of the U 2 M 0 compound in the Cl U structure [18]. In addition, there is one theoretical attempt 
to model some structural, elastic and thermal properties of this compound using using density-functional theory (DET) 
and the quasi-harmonic Debye theory [12]. 

It is important to note here that all previous works assumed the Cl U structure as the ground state of the U 2 M 0 
compound. The results of this work show that this compound is unstable under the 05 deformation related to the C 66 
elastic constant. Consequently, the use of an unstable compound in the ECI expansion to describe the ground state of 
the U-Mo system could be the source of non-physical conclusions. This result reveals the importance of an accurate 
characterization of the physical and thermal properties of this compound and, mainly, the ground state of a system 
previously to apply the cited theoretical methodologies. 

The purpose of this work is to study the structural and electronic properties of U 2 M 0 to search for the origin of 
the instability under the £>6 deformation. The paper is organized as follows. Section 2 gives the methodology used to 
compute the elastic constant. Section 3 presents the analysis of the structural instability observed in U 2 M 0 together 
with the electronic properties as DOS, bands and Eermi surface. 

2. THEORETICAL METHODOLOGY 

2.7. Structure optimization and calculation parameters 

The assumed crystal structure for the ground state of the U 2 M 0 compound is the body-centred tetragonal struc¬ 
ture with the space group 14 /mmm (#139), MoSi 2 type prototype with Strukturbericht designation Cnh [19]. The 
molybdenum atoms are at the positions 2(a ): (0,0,0) and (V 2 , V 2 , V 2 ), and the uranium atoms at the positions 4(e ): 
(0,0, z), (0,0, -z), (V 2 , V 2 , V 2 + z) and (V 2 , V 2 , V 2 - z) ; where z = 0.328 + 0.002. The experimental lattice parameters 
are ao = bo = 3.427A, cq = 9.834A [20]. 

The total energy of U 2 M 0 was computed using the full potential LAPW method based on density functional 
theory as it is implemented in the WIEN2k code [16]. This code uses the full-potential LAPW-rlo method that makes 
no shape approximation to the potential or density. The electronic exchange-correlation interactions were treated 
within the generalized gradient approximation of Perdew, Burke and Ernzerhof [22], as no experimental or theoretical 
evidence of strong correlations was found on the system studied here. The radii of the atomic spheres (Rmt) selected 
for U and Mo were RmTu = ^.4 a.u. and RmTmo - order to describe the electronic structure of all the atoms 

and its orbitals, in this calculation, the APW + lo basis set was selected. This choice was made considering the fact 
this method allows to reach converged values at a lower Rkmax, and thus, the corresponding computational time is 
drastically reduced. Local orbital extensions were included to describe the semicore states by means of the method 
(APW + lo ), as it is implemented in the WIEN2k program. The cutoff parameter that controls the convergence in the 
expansion of the solution to the Kohn-Sham equations is chosen to be Rkmax = 8 , where kmax is the plane wave 
cutoff and RMT is the smallest of all the atomic sphere radii. The maximum I values for partial waves used inside the 
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Figure 1: The crystal structure of the assumed ground state for the U 2 M 0 compound [19, 21]. 



Figure 2: Comparison of the density of states around Ep with and without spin orbit. k-points= 10000. 


atomic spheres and for the non-muffin-tin matrix elements were selected to be Imax =10 and Imax = 4, respectively. 
The charge density cut-off Gmax was selected as A mesh of 10000 k-points was taken in the whole Brillouin 

zone, k-space integration is done by the modified tetrahedron-method [23]. The iteration process is repeated until the 
calculated total energy converges to less than \x\Q~^Ryjcell, and the calculated total charge converges to less than 
1x10 - 6 elcelL A complete determination of the tetragonal structure requires the determination of three parameters 
(a, b and the internal position parameter). In order to calculate the internal parameters of the crystal structure we 
employed the mini LAPW script implemented in the Wien2K package . 

The theoretical lattice parameters at 0°K obtained were a = b = 3.4437A, c = 9.6842A and z = 0.3216, in 
good agreement with experimental results. The calculations were performed without including spin-orbit interactions 
because no significant effects were observed near the Fermi level, as it is shown in Fig. 2 for the total DOS. 

2.2. Elastic constants C 44 and Cee calculations in U 2 M 0 

The elastic constants are obtained by computing the second order derivatives of the total energy with respect to the 
distortion parameter. The elastic constants of U 2 M 0 were studied in a previous work [12], using the density functional 
theory implemented in the WIEN2k code. These calculations were performed by using a package included in the code 
which calculates elastic tensors of tetragonal phases [24]. The results are presented in Table 1. 

The elastic constants are very sensitive to the data point selection due to numerical inaccuracies of the calculations 
around the minimum of the total energy. In order to have a better control of the outcome in the cases of constants with 
very small values, such as C 44 and C 66 , the calculations of the total energies were repeated using the Wien2k code 
but not its elastic packages. This was achieved by considering the structure in the primitive triclinic unit cell with 
space group PI (see Fig. 3 ). After global relaxation, the primitive lattice vector sive a = b = c = 5.4199A and angles 
a = JS = 142.954° and y = 53.395°. The basis comprises only three atoms, one Molybdenum at position (000), one 
Uranium at (u u 0) and the other Uranium at (u u 0); being u = 0.3216. All the positions are referred to the primitive 
lattice vectors a, b and c in crystal parameters. See Fig. 3. 
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Table 1: Calculated elastic constants (in GPa), bulk modulus B, Young’s modulus E, shear modulus G and Poisson’s 
ratio y. 


Cii 

254 

B 

180 GPa 

Gi2 

161 

E 

102 GPa 

Gi3 

125 

G 

36 Gpa 

G 33 

295 

y 

0.4 

C 44 

38 



^66 

~20 





Figure 3: U 2 M 0 primitive unit cell. It is a triclinic crystal structure with vector modules: |a| = |b| = |c| = 5.4199A and 
angles: a=p = 142.953766°, y = 53.394693°. 


The elastic constants C 44 and are calculated through the volume conserving monoclinic shear distortions £>4 
and £> 6 , respectively. Matrices for each distortion are given by Eq. 1 and 2, 
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where, 6 is the dimensionless distortion parameter. 

The energy associated to this distortion can be written in the following way: 

E{V, 5) = E(Vo, 0) + yo(2T45 + 2 c44<5^) , (3) 


E(V, 6) = EiVo, 0) + VoilTeS + IceeS^), (4) 

where Vo is the volume of the unstrained system, £’(Vo, d) is the total energy computed by ab initio codes in this work, 
T 4 and T 6 are the elements of the stress tensor in Voigt notation (T 4 = re = Txy). Therefore, the elastic constants 
are deduced from the following expressions: 
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3. RESULTS AND DISCUSSIONS 


3.1. Structural instability in U 2 M 0 



Figure 4: Total energy for different values of the distortion parameter d in U 2 M 0 after the relaxation of internal 
positions. Spin-orbit coupling was not included in the calculations. 

The low magnitude of the C 44 and constants led the question whether a soft phonon mode is related to those 
deformations or if a more precise calculation around the minimum of energy could reveal hidden structural anomalies 
due to strong electron-phonon coupling. 

Following the procedure outlined in section 2, a lattice instability was found for the £>5 deformation but none for the 
D/[ one. Fig. 4 shows the total energy versus the distortion 6 for the compound U 2 M 0 in CIU structure. The total 
energy reaches a minimum value for a displacement of d = +0.03. The magnitude and shape of the total energy vs. 
distortion curve obtained under the application of the £>6 distortion in U 2 M 0 is very similar to the one calculated in 
the a-U system, where a charge density wave is observed [25]. However, whereas the CDW in a-U is spontaneous, 
the effect observed in U 2 M 0 is induced by deformation and could be a remnant of the effect observed in pure U. If 
so, U 2 M 0 will be another example of a three-dimensional system with this kind of phenomenology. In addition, the 
instability remarks that the ground state of U 2 M 0 compound is a crystalline structure different to the Cl U. 

3.2. £>6 distortion in U 2 X (X= Tc, Rh, Ru, Cr, Pd, W, Al, Nb, Z) metastable compounds 



(a) X element in the same period (5). 

Figure 5: Comparison of the total energy variation under 
U 2 X compounds. Eq corresponds to the value of energy at 



(b) X element in the same group (VI). 

application of the 06 distortion parameter 6 for different 
the undistorted lattice. 


In order to study the effect produced by the occupation of the electronic states on the structural instability related 
to the £>6 distortion, the behaviour of the original U 2 M 0 structure was analysed by replacing the Mo atoms with U 
atoms (U 2 U) and also with other nearby elements of the periodic table (Tc, Rh, Ru, Cr, Pd, W, Al, Nb, Zr). The 
calculations will give information whether it is a specific effect related to only the U 2 M 0 compound or is a more 
general effect involving other elements with different occupation of the d-states. The £>6 distortion matrix was applied 
to the unit cell (see Fig. 3) after volume optimization and internal parameters relaxation of each structure in the space 
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group MImmm. The total energy differences relative to the values obtained for the undistorted cells as a function of 
the distortion parameter S, are shown in Fig. 5a and 5b. The former compares elements situated in the same period 
(5) and the later one compares elements on the same group (VI). 

Whereas it is not observed any instability in the U 2 U compound, the existence of such instability is observed in 
other U 2 X compounds. The result emphasizes two facts: i) the U atoms involved in the CDW in a-U behave different 
in the U 2 X compounds and, ii) the interaction and number of d-electrons play a crucial role in the existence of the 
instability. Figs. 5a and 5b show that the structural instability is observed only if the transition elements have at least 
five d-electrons. The results reveal that the interaction between d- and f-electrons plays a fundamental role in the 
instability. Therefore, the effects of the hybridization should be observed in the electronic properties such as DOS, 
band and Fermi surface. 

3.3. Electronic properties analysis 
3.3.1. DOS andprojected-DOS 

The first electronic properties computed were the DOS and partial DOS. These calculations were done relativis- 
tically, without including spin-orbit interaction, as spin-orbit does not seem to have a large effect, as it was shown in 
Fig. 2. A comparison of the partial DOS is given in Fig. 6a for the undistorted (d=0) and in Fig. 6b for the distorted 
(d=0.03) cases. A small change of the projected f-DOS is observed below and around the Fermi level. It is clear from 
the projected DOS that the f-states of U play a fundamental role near the Fermi level (FL), similarly to a-Uranium. 




Figure 6: Electronic density of states and projections onto the LAPW spheres. k-points=50000 


3.3.2. Electronic energy band structure 

The electronic band structure for U 2 M 0 was calculated through the path shown in Fig. 7. The orbital character is 
illustrated in Figs. 8a and 8b for the undistorted and distorted structure, respectively. The f-states of U and the d-states 
of Mo are plotted in a 1 eV range around EF. 



Figure 7: Brillouin zone of ORCE 2 (Face-centred orthorhombic 1 /a^ < 1/b^ l/c^) type lattice [26] with the selected 
path for bands calculations. 

The bands plot for the undistorted structure shows the Fermi level (FL) intersected mainly by the f-states of U. 
However, there are zones near the FL showing d-states of Mo atoms hybridized with the f-states. They are located 
near the X, Y and Z points of the Brillouin zone. They are coloured in purple in Fig. 8. For the distorted case, it is 
observed in Fig. 8b an important rearrangement of the bands. The lowered symmetry of the distorted structure breaks 


6 
























the degeneracy of the f-bands of U. There is an important split of the bands just below the FL near the F point. The 
main splits are observed in the directions F-X, F-Y, F-Z and F-L. Those changes modify completely the number of 
states crossing the FL modifying the structure of the Fermi surface. It is still observed the hybridization of d-Mo and 
f-U states, but those states move below the Fermi level near the Y, X and Z points. Consequently, there are occupied 
states in the distorted structure with a lower energy than the same states of the symmetric structure. Therefore, the 
low symmetry structure is stabilized by the relocation below the FL of the hybridized states between Mo and U atoms 
together with the f-band split of U, as shown in Fig. 8b. At variance with a-U, the d-sates of Molybdenum play an 
important role in stabilizing the distorted structure. Changes like those observed in the bands of the distorted U 2 M 0 , 
stabilizing structures with lower symmetry, are usually referred to as a Peierls distortion. However, whereas the Peierls 
distortion in a-U is spontaneous, in U 2 M 0 is induced by the deformation D^. 

The Cllb structure can be interpreted as formed by parallel linear chains along the z-directions each one composed 
by U-Mo-U blocks. The changes due to the distortion modify the distance between the Mo and U atoms. The 
hybridization mainly due to interaction inside the U-Mo-U blocks is slightly modified. Those states just move below 
the Fermi level contributing to stabilize the distorted structure. The changes in distance among chains modify the U-U 
interaction and reproduce the same situation observed in a-U, i.e. the split of f-state. The induced Peierls distortion 
in U 2 M 0 arises from this mechanism. 



Figure 8: Band structure (left) and predominant character acting nearby the Fermi level (right) obtained following the 
path: Y-r-X-L-r- Zk-points=50000. 


3.3.3. Fermi surface 

The Fermi surfaces for the undistorted and distorted systems are shown in Figs. 9 and 10 respectively. There is 
a clear topological change in the Fermi surface close to the X point, and mainly, in the F-X, F-Y and F-Z directions. 
The band # 30 has disappeared in the distorted structure. In addition, the band #31 becomes disconnected along the 
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(a) Band # 30 (b) Band #31 (c) Band # 32 (d) Band # 33 

Figure 9: Fermi surfaces corresponding to the different intersecting bands with d = 0. k-points=50000. 



(b) Band # 32 (c) Band # 33 

Figure 10: Fermi surfaces corresponding to the different intersecting bands with 6 = -0.03. k-points=50000. 




F-X direction. The states along the other two directions produce a connected Fermi surface, observed clearly in band 
#33. The mechanism is similar to the Peierls distortion that occurs in systems with narrow bands intersecting the 
Fermi level, as was explained by Soderlind for a-\] [27]. Both conditions, necessary for the existence of a charge 
density wave, are met in the case of the structure of the compound Cl 1^. 

4. Conclusions 

Besides the technological importance of the U-Mo solid solution in nuclear energy research, there are few theo¬ 
retical attempts to model this system. All of previous works assumed the Cl U structure as the ground state of the 
U 2 M 0 compound. This work reports the existence of a structural instability of this compound under the deformation 
06, related to the elastic constant. This result will be helpful in the modelling of the U-Mo system with theoretical 
methods using ground state information. The structural and electronic properties of U 2 M 0 are studied in this work to 
find the source of the instability under the 06 deformation. It is found that the split of the f-states of U together with 
the energy lowering of the hybridized states, mainly between d-Mo and f-U states, are responsible of the stabilization 
of the distorted structure. The Fermi surface is analysed. It is found a clear topological change in the Fermi surface 
close to the X point, and mainly, in the F-X, F-Y and F-Z directions. Changes like those observed in the bands of the 
distorted U 2 M 0 , i.e. the split of f-bands and the lowering of f-states, stabilizing structures with lower symmetry, are 
usually referred to as a Peierls distortion. However, whereas the Peierls distortion in a-\] is spontaneous, in U 2 M 0 is 
induced by the deformation It remains to study whether a connection between the £>6 deformation and a CDW 
can be established and to find the structure of the ground state of the U 2 M 0 compound. 
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